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SYMBOLS 


area  of  tube 

th 

d.  mean  diameter  of  i  tube 

fK 

critical  buckling  load  for  i  "  tube 
maximum  wing  tip  force 

F^^  that  portion  of  F^,  required  to  rotate  wing  quasi-statically, 

i.e.  ,  required  to  collapse  wing  connection  structure 

F^^  that  portion  of  F^  required  to  impart  angular  momentum 

to  wing  mass 

moment  of  inertia  of  i^^  tube 
mass  moment-of-inertia  of  wings 

K  elastic  modulus  of  receiving  medium  (force /unit  length 

compression) 

Lj  length  of  i^^  tube 

m  total  mass  of  aircraft 

m'  mass  of  aircraft  excluding  engine 

m^^  mass  of  occupant 

m^  mass  of  wings 

Pj  longitudinal  force  required  to  collapse  first  bay 

r  distance  from  center  of  rotation  to  the  point  of  application 

8  length  of  first  bay 

s  total  displacement  of  center  of  mass  of  aircraft  less 

engine  during  impact  (s'  =  x  +  s) 

t  duration  of  impact 


IV 


u 


g 


V 


V 


1 


V 


o 


X 


Xi.  Vi,  Zi 
xi>  yi.  Zi 


Y 

z 

A 

Axi 


tube  wall  thickness,  tube 

energy  required  to  collapse  (completely)  the  main  gear 
initial  velocity  of  aircraft 

velocity  of  fuselage  upon  its  contact  with  ground 
initial  vertical  velocity  of  aircraft 
maximum  deflection  of  receiving  medium 
coordinates  of  forward  end  of  i^^  tube 
coordinates  of  aft  end  of  i^^  tube 
a  dimensionless  variable 
a  dimensionless  variable 

reduction  in  length  of  the  i^^  tube  (associated  with  the 
assumed  mechanism) 

total  deflection  of  seat>cusMon  combination  associated 
with  a  critical  vertical  deceleration  of  occupant 

angle  of  wing  rotation 

direction  cosine  of  i^^  tube  for  the  longitudinal  direction 
yield  stress 

maximum  angixlar  velocity  of  wings 
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SUMMARY 


This  report  covers  the  development  of  limited  mathematical  models  to 
predict  in  qualitative  terms  (see  Limitations  of  Analysis,  page  24)  the 
crashworthiness  of  a  light  airframe  structure  (of  the  Piper  J-3  type) 
under  three  impact  conditions. 

Two  indices  of  crashworthiness  were  selected:  (1)  the  aircraft  velocity 
required  to  cause  structural  collapse  short  of  general  cabin  collapse, 
and  (2)  the  maximum  deceleration  applied  to  the  occupants  during  the 
crash. 

Engineering  analysis  and  physical  principles  are  applied  to  each  appro¬ 
priate  mathematical  model  to  develop  the  crashworthiness  indices. 
Sample  calculations  were  made  for  the  Piper  J-3  in  three  typical  crash 
impact  configurations:  nose  impact,  nose  and  wing  impact,  and  impact 
on  the  landing  gear  (considering  only  the  vertical  velocity  component). 

The  preliminary  study  reveals  that  the  crashworthiness  indices  selected 
are  potentially  useful  in  evaluating  structural  modifications  such  as  the 
eight  items  discussed. 
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CONCLUSIONS 


1.  This  preliminary  study  indicates  that  the  crashworthiness  indices 
selected  are  potentially  useful  in  evaluating  structural  modifications 
to  increase  the  crashworthiness  of  a  given  aircraft  structure. 

2.  The  efficiency  of  energy  absorption  of  the  first  bay  has  a  significant 
effect  upon  the  selected  indices  and  hence  upon  crashworthiness. 

3.  Considerable  promise  of  improved  crashworthiness  appears  probable 
from  changes  in  cabin  structure.  If  the  cabin  were  made  much 
stronger  than  other  airframe  structure,  impact  forces  would  be 
transmitted  through  the  cabin  and  the  energy  absorbed  in  deforma¬ 
tion  elsewhere. 

4.  This  preliminary  study  indicates  that  a  more  extensive  investigation 
of  the  collapse  mechanism  involving  computer  simulation  and  time 
sequence  is  feasible. 
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RECOMMENDATIONS 


It  is  recommended  that: 

1.  Suitable  experiments  be  designed  and  conducted  to  check  the 
validity  of  the  theoretical  assumptions,  thus  removing  some 
of  the  current  information  gaps, 

2.  For  a  more  general  treatment,  three  approaches  be  explored 
in  order  to  determine  which  best  serves  the  study  objectives; 
upon  selection  of  the  best  approach,  it  be  adopted  for  a  general 
computer  program.  The  three  approaches  are  described  as 
follows; 

a.  Separate  analyses  would  be  developed  for  aircraft  of 
basically  different  structures  with  separate  subroutines 
for  different  substructures.  If  a  final  computer  program 
designed  to  analyze  any  light  aircraft  were  to  evolve  from 
this  approach,  logical  branching  to  alternate  routines  and 
subroutines  would  be  required. 

b.  A  second  approach  to  generality  would  entail  a  basic  de¬ 
scription  of  any  airframe  element  (a  tube,  a  monocoque 
bay,  etc.)  in  terms  of  fundamental  properties  (geometric 
and  structural  behavior  properties).  Then,  by  locating 
these  elements  with  respect  to  one  another,  euiy  airframe 
could  be  mathematically  simulated.  A  simulated  impact 
configuration  would  be  applied  to  such  a  structure  on  the 
computer.  Then,  employing  the  physical  principles  in¬ 
volved,  the  collapse  mode  and  force  transmission  would 
be  determined, 

c.  A  third  approach,  of  possible  usefulness,  would  involve  an 
abstraction  of  structural  behavior  from  the  external  impact 
conditions.  Based  upon  statistical  accident  data  and  engi¬ 
neering  analysis,  a  number  of  potential  collapse  modes  or 
mechanisms  would  be  selected  for  examination.  For  each 
of  these  mechanisms,  crashworthiness  indices  would  be 
developed,  such  as  specific  energy  absorption  or  a  specific 
Ci'itical  force  index.  A  meaningful  computer  program  that 
followed  this  approach  would  compute  the  crashworthiness 
indices  for  a  given  structure  and  would  also  provide  statis¬ 
tical  norms  and  ’’significance"  weight  factors  for  each  of 
the  indices  computed.  The  index  norm  and  weight  factor 
would  conceivably  depend  upon  the  design  mission  and  per¬ 
formance  features  of  the  particular  subclass  of  aircraft. 
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STATEMENT  OF  PURPOSE 


The  objective  of  the  current  study  is  to  determine  significant  indices 
of  crashworthiness  as  related  to  the  structure  oi  light  aircraft.  The 
question  for  which  an  answer  is  sought  is  simply:  How  may  proposed 
changes  in  a  given  structure,  intended  to  increase  crashworthiness, 
be  evaluated  most  meaningfully?  It  is  not  intended  that  these  indices 
be  used  to  rate  crashworthiness  of  various  makes  and  models  of  air¬ 
craft  relative  to  one  another  (as  other  factors  would  enter  into  such 
a  rating)  but  rather  that  these  indices  should  serve  to  point  directions 
of  potential  improvement  in  a  given  structural  design.  Perhaps  useful 
generalizations,  such  as  a  recognition  of  the  critical  nature  of  energy¬ 
absorbing  properties  of  the  first  bay,  will  emerge  as  the  study  pro¬ 
gresses. 
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PHILOSOPHY  OF  ANALYSIS 


Survivability  or  nonsurvivability  of  a  lightplane  accident  is  generally 
determined  by  the  degree  to  which  the  aircraft  structure  around  the 
occupants  collapses  or  retains  its  shape.  In  an  accident  in  which  the 
cabin  structure  collapses  to  a  fraction  of  its  normal  volume,  the  chances 
of  occupant  survival  are  correspondingly  reduced.  Additionally,  if 
decelerations  are  great  enough  to  exceed  the  limits  of  the  occupant 
restraint  systems,  severe  injury  could  also  result. 

The  preliminary  study  employs  as  crashworthiness  indices:  (1)  the 
maximum  impact  velocities  which  can  be  sustained  by  the  aircraft  in 
various  impact  configurations  without  causing  general  cabin  collapse, 
or  in  the  case  of  vertical  motion,  without  causing  spinal  injury;  and 
(2)  the  decelerations  associated  with  these  limiting  velocities. 

Nearly  all  light  aircraft  in  use  today  are  of  two  types  of  construction: 

(1)  welded  steel  tubing  with  fabric  covering,  or  (2)  all>metal  semi- 
monocoque.  The  first  portion  of  this  study  deals  with  tube  and  fabric 
structure,  using  the  Piper  J-3  as  an  example. 

As  it  is  not  feasible  to  consider  all  possible  accident  impact  config¬ 
urations,  a  representative  few  are  selected  for  analysis  that  typify 
conditions  most  frequently  found  in  actual  accidents.  A  survey  of  a 
representative  group  of  survivable  Piper  J-3  accidents  gives  the 
following  distribution  of  impact  configurations  classified  on  the  basis 
of  the  contact  points  during  the  primary  impact. 


IMPACT  CONFIGURATIONS  OF  116  PIPER  J-3  ACCIDENTS 


No.  of 

Percent 

Impact  Configuration 

Cases 

of  Cases 

Nose  only 

33 

28.45 

Nose  and  wing 

31 

26.  72 

Nose  and  main  landing  gear 

22 

18.97 

Nose,  landing  gear  and  wing 

13 

11.21 

Landing  gear  only 

8 

6.89 

I.anding  gear  and  wing 

5 

4.31 

Wing  only 

4 

3.45 

116 

100.00 
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Using  the  foregoing  table  as  a  guide,  mathematical  models  of  impact 
configurations  are  employed  with  the  same  primary  contact  points.  It 
is  assumed  that  the  flight  path  is  along  the  longitudinal  axis  of  the  air¬ 
craft,  except  the  "gear  only"  configuranon.  An  ideal  elastic  receivi. 
medium,  simulating  a  typical  soil,  is  assumed  for  purposes  of  energy 
absorption  calculations. 


ANALYTICAL  METHOD 


To  estimate  the  critical  velocity  for  a  given  impact  configuration, 
that  is,  the  velocity  of  the  aircraft  required  to  cause  structural  collaps 
short  of  general  cabin  collapse,  an  energy  approach  is  employed.  The 
kinetic  energy  of  the  aircraft  just  prior  to  Impact  is  equated  to  the 
energy  dissipated  in  various  forms  during  the  impact.  This  latter 
includes  the  energy  to  compress  the  receiving  medium,  the  work  of 
plastic  collapse  of  aircraft  structure,  and  the  energy  dissipated  in 
internal  damping  throughout  the  aircraft  (other  than  that  of  visible 
plastic  deformation). 

To  determine  the  maximum  deceleration  experienced  by  the  occupants, 
the  second  index  of  crashworthiness,  a  force  transmission  analysis  is 
considered  in  which  the  forces  required  to  cause  each  stage  of  the 
progressive  collapse  arc,  computed.  The  maximum  force  .found  may  be 
used  to  obtain  the  maxim un>  deceleration  experienced  if  the  effective 
mass  of  the  aircraft  is  known. 

If  the  mechanism  or  mode  of  collapse  is  not  evident,  it  may  be  useful 
to  bracket  the  unknown  collapse  force  by  upper-bound  and  lower-bouJtid 
values.  In  particular,  for  an  assumed  mechanism,  the  load  required 
to  sustain  the  mechanism  will  be  greater  than  or  equal  to  the  actual 
collapse  forces;  hence,  the  load  will  form  an  upper  bound.  Corre¬ 
spondingly,  for  an  assumed  stable  internal  force  distribution  satisfying 
equilibrium,  the  associated  loads  are  always  less  than  or  equal  to  the 
required  collapse  forces,  constituting  a  lower  bound. 

/s  a  practical  computational  technique,  a  probable  collapse  mode  is 
first  assumed  which  would  yield  an  upper  bound.  This  then  would  be 
examined  from  a  force  distribution  viewpoint  to  determine  whether  it 
also  satisfies  the  lower-bound  requirements.  If  not,  a  lower  bound  can 
be  obtained  simply  by  a  prop>ortionate  reduction  in  force  magnitudes  to 
the  i>oint  where  the  lower-bound  requirements  are  satisfied.  If  the 
resulting  upper  and  lower  bounds  are  sufficiently  close  to  each  other, 
their  average  constitutes  a  satisfactory  approximation  to  actual  collapse 
forces . 
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APPLICATION  OF  METHOD 


Three  impact  configurations  have  been  considered  initially  for  which 
crashworthiness  indices  are  developed  in  terms  of  structural  para¬ 
meters  (tube  and  fabric  structure  only,  to  date).  These  three  con¬ 
figurations  are:  (1)  nose -only  contact,  (2)  nose-and-wing  contact, 

(3)  gear -only  contact,  considering  vertical  velocity.  The  detailed 
analysis  of  each  of  these  is  treated  separately  as  follows. 

NOSE -ONLY  IMPACT  CONFIGURATION 

A  simplified  mathematical  model  using  lumped  parameters  is  present¬ 
ed  schematically  in  the  following  diagram: 


RECEIVING 

MEDIUM 

(ELASTIC 


MASS  BAY 
(RIGID)  (COLLAP¬ 
SIBLE) 


AIRCRAFT 
LESS  ENGINE 


SPRING) 


The  pattern  of  impact  and  collapse  is  assumed  to  entail  a  compression 
of  the  receiving  medium  initially,  with  the  reaction  force  increasing 
linearly  to  the  point  where  the  critical  axial  load  causes  a  complete 
crushing  of  the  first  bay  (provided  the  forces  required  are  less  than 
the  critical  forces  for  cabin  collapse). 


Symbol  Notation: 


elastic  modulus  of  receiving  medium  (force /unit  length 
compression) 


V  initial  velocity  of  aircraft 


X  maximum  deflection  of  receiving  medium 
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P|  longitudinal  force  required  to  collapse  first  bay 

s  length  of  first  bay 

m'  :  lass  of  aircraft  excluding  engine 

m  total  mass  of  aircraft 


Considering  the  entire  aircraft  ir.  an  energy  balance,  the  kinetic  energy 
prior  to  impact  would  equal  the  energy  required  to  compress  the  re¬ 
ceiving  medium,  the  energy  required  to  crush  the  first  bay,  and  the 
energy  dissipated  in  internal  damping  throughout  the  aircraft-  The 
latter  quantity  is  of  a  nebulous  nature,  although  dependent  in  part  upon 
the  magnitude  of  Pi  and  the  length  s.  As  an  operating  expedient,  the 
average  axial  force  transmitted  through  the  first  bay  is  taken  as  half 
of  P  (the  critical  first-bay  load  under  static  conditions),  although  it 
is  known  that  a  column  force  falls  off  markedly  after  buckling  has 
started.  This  simplifying  assumption  of  treating  the  average  force  as 
half  the  maximum  static  force  becomes  plausible  when  one  considers 
the  nebulous  nature  of  internal  damping  energy,  the  existence  of  ob¬ 
structions  to  the  free  column  buckling  of  the  tubular  elements,  the 
direct  axial  resistance  offered  by  nonstructural  elements  (sheet-metal 
covering,  fuel  tank,  etc.),  and  the  dynamic  column  load  over-shoot 
associated  with  high  deformation  rates.  Thus  the  energy  balance  may 
be  written: 


2 

/ 

initial 

kinetic 

energy 


\ 

X 

receiving  aircraft 
medium  deformation 

energy  energy 


(1) 


In  a  similar  fashion,  an  energy  balance  may  be  written  for  the  aircraft 
without  engine,  in  terms  of  motion  of  its  mass  center.  The  principle 
used  may  be  stated  as  follows;  The  work  of  external  forces  displaced 
through  the  mass  center  displacement  equals  the  change  in  kinetic 
energy  associated  with  the  velocity  of  the  mass  center.  Thus,  for  an 
average  external  force  and  a  total  displacement  of  mass  center 
equal  to  (x  +  s),  2 
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-  Y  (x  s) 


2 


p  change  in 

negative  work  of  —  kinetic 

energy 


(2) 


Solving  equation  (2)  for  x  and  substituting  into  (i)  yields 


This  may  be  simplified  to 


Z2-  (2  ♦  ^  Sj.)z  ♦  (1  ♦  ^)  *  0  (3) 

where  Z  is  a  symbol  for  the  dimensionless  qviantity 

Z  .am.  '-‘5 

Pi* 

If  Pi  is  first  determined,  then  by  means  of  equations  (3)  and  (4),  the 
critical  velocity  v  may  be  computed. 

For  the  computation  of  Pi ,  the  longitudinal  force  required  to  collapse 


the  first  bay. 

the  following  notation  is  introduced; 

^i»  Yi*  *i 

coordinates  of  forward  end  of  i^^  tube 

M,  yi,  zi 

coordinates  of  aft  end  of  i*^  tube 

H 

tube  wall  thickness,  i^  tube 

di 

mean  diameter  of  i^^  tube 

Li»  Aj,  li 

length,  area,  and  moment  of  inertia  of  i^^  tube 

^cri 

critical  buckling  load  for  tube 

Xxi 

direction  cosine  of  i^^  tube  for  the  longitudinal 
direction 

•y 

yield  stress 

I 


For  a  sufficiently  slender  pinned-end  column,  the  critical  load  is 
given  by 

»  »^EI 

p  m  * 

cr  j2  * 

and  for  a  fully  fixed-end  colximn, 

F  -iiis.. 

cr  ^2 

In  a  welded  tubular  structure,  some  intermediate  end  fixity  condition 
exists;  hence, 

2s2bi, 


is  taken  as  a  plausible  value  for  the  buckling  load  of  a  slender  coltunn. 

A  column  of  low  slenderness  ratio,  on  the  other  hand,  would  undergo 
plastic  rather  than  elastic  buckling,  the  critical  load  for  which  may  be 
crudely  approximated  by  yield  stress  times  the  cross-sectional  area. 


For  circular  steel  tubes. 


Fcr^  *  OyAi  • 


E  ■  30  X  10^  ksi 


I  * 


A  »  w(d)t  . 


Equations  (5)  and  (6)  thus  reduce  to 

''crj  •  23-2 

Fcr^  »  Oyirtidi 

where  the  smaller  value  of  F  would  govern. 


From  geometry 


k  ■  Ji^i  -  * 

The  longitudinal  component  of  F^.^.  is  obtained  by  multiplying  by 
direction  cosine,  ^vhere  X^i  is  given  by 

Xj  - 

^xi  * 

Consequently,  the  total  critical  longitudinal  load.  Pi,  is 

Pi  “  I  Pi^xi  •  (10) 

An  example  of  computations  for  the  nose -only  type  of  impact  is  given 
in  the  Appendix,  pages  27,  28,  and  29. 

NOSE-AND-WLNG  CONTACT  CONFIGURATION 

A  schematic  mathematical  model  of  the  aircraft  for  this  impact  con¬ 
figuration  is  presented  below. 


AGENT 
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The  symbol  notation  employed  for  the  analysis  of  this  configuration 
appears  below  (previously  introduced  notation  would  still  apply). 

maximum  wing-tip  force 

that  portion  of  F^  required  to  rotate  wing  quasi-statically, 
i.e.  ,  required  to  collapse  wing  connection  structure 

F^2  portion  of  F^  required  to  impart  angular  momentum 

to  wing  mass 

m^  mass  of  wings 

mass  moment-of-inertia  of  wings 

r  distance  from  center  of  fuselage  to  wing  tip  (point  of 

force  application) 

w  maximum  angular  velocity  of  wings 

"S'  total  displacement  of  center  of  mass  of  aircraft  less 

engine  during  impact  (s'  =  x  +  s) 

t  duration  of  impact 

Y  a  dimensionless  variable 

Simplifying  assumptions  made  in  the  ensuing  analysis  are  as  follows; 

a.  Under  the  crash  kinematics,  the  wing  acquires  an  angular 
momentum  where  the  maximum  angular  velocity,  «  ,  is 
approximated  by 


V  still  denoting  the  velocity  of  the  aircraft  prior  to  impact. 

b.  The  average  wing  force  during  the  impact  time  interval  is 
assumed  to  be  half  of  the  maximum  force,  F^  . 

c.  The  average  velocity  ^ne  aircraft  center  of  gravity 
(excluding  engine)  during  the  impact  is  assumed  to  be  half 
of  the  initial  velocity,  v. 
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d.  The  fuselage  does  not  rotate  due  to  eccentric  loading  of 
the  wing,  but  continues  its  free-path  movement  linearly. 

e.  The  wing  contact  agent  is  fully  rigid  and  therefore  does 
not  deflect. 

Considering  an  energy  balance  associated  with  motion  of  the  center 
of  mass  of  the  aircraft  less  engine,  we  have 

8*4  a*v^  . 

2 

Then,  since  ~  +  F^,2» 


^  ^V2  _  1,2 

*  *  2  “  ' 


However,  s 


work  relative  to  the  center  of  mass  of  the  wing 
of  force  Fw2  upon  the  wing; 
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From  kinematics,  the  time  may  be  found: 


8  *  j  t 


t  * 


2«»v 


Pi  +  P, 


vl 


(12) 


The  angular  impulse /angular  momentum  relationship  may  be  used  to 
determine  the  force  ' 

Angular  Impulse  *  rPAt  *  rPt 

t  *  At  *  (tj*  -  t^) 

Iq  ■  time  initial 
tf  ■  time  final 


Angular  Momentum  =  Iw. 

Since  the  angular  impulse  of  ‘he  external  forces  acting  on  a  system  of 
particles  with  respect  to  an  ;xxis  fixed  in  space  during  any  time  interval 
is  equal  to  the  change  of  angular  momentum  of  the  system  with  respect 
to  the  same  axis,  the  following  expression  may  be  written: 


Frt  ■  I  ( uif  «  ) 

«  U 

■  0 

Frt  «  I<D  . 


f’v2  ^  V 

tr  »  -  "f”  r 


F  «  ^ 
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Thus,  substituting  equation  (12)  and  simplifying, 

P  ♦  P 

«  _  Wi 

N 

An  energy  balance  for  the  entire  aircraft  may  be  employed  to  determine 
the  critical  initial  velocity,  v: 

sd  ,  +  V  ^  Lai  ^  idl  (14) 

2  2  2  2  2 

initial  work  on  energy  of  energy  of  kinetic 

kinetic  receiving  deformation  deformation  energy 

in  first  bay  for  wing  rotation  of  wing 


energy 


medium 


Substituting  for  x,  and  w  the  expressions 


X  «  I  -  s  , 


■w^ 

■*v  3 


equation  (14)  becomes 


(1  -  ~)  ■  K  (s^  -  288  t^)  ♦  P|8  F^ji  . 

Noting  that  ji  “  j  *  jl  -  ^  j  ,  equation  (11)  may  be  introduced 


into  the  left-hand  side  of  the  above: 


(?!  ♦  Fvi)  i  ■  K  (i^  -  2i8  s^)  ♦  P^s  ♦ 


-  (2  ♦  5- -  Isa]  Y  ■>  In-  ill.  0 

1  ■'  K8  K8  j  I  Ks  J 
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where 


Y 


■fv2 


(16) 


Equation  (15)  can  be  solved  for  the  dimensionless  quantity  Y  (assuming 
Pi  and  Fwi  are  known  quantities  at  this  point).  The  critical  velocity  v 
can  be  obtained  then  from  expression  (16)  for  Y. 

Turning  attention  next  to  the  computation  of  (the  Pi  computation 
having  previously  been  indicated,  expression  10),  a  plausible  mechanism 
of  airframe  collapse  at  the  wing  attachment  may  be  considered.  One 
such  mechanism  is  a  rigid-body  rotation  (relative  to  the  fuselage)  of 
the  wings  and  top  of  the  cabin  (between  the  wings)  about  the  geometric 
center  of  the  attachment  points.  Under  such  a  rotation,  deformation 
occurs  in  the  fuselage  bays  adjacent  to  the  wing  attachments,  the 
pattern  of  which  is  determined  from  geometry  and  from  the  assxunption 
of  no  extension  of  any  member.  The  energy  of  deformation  is  then 
computed  as 


i 


energy  of  ’^formation 


where 

F^r^  critical  column  lead  for  the  i^^  tube 

1 1  is  the  reduction  in  length  of  the  i^^  tube  (associated  with  the 
assumed  mechanism) 

Equating  the  external  work  done  by  the  applied  force  Fwi  to  the 
energy  of  deformation,  we  obtain 

Fyjr60  -  IPcrj^i 

or 

I 

*^.1  ■  — TiT"  (17) 


where 

r  is  the  distance  from  the  center  of  rotation  to  the  point  of 

application  of  Fvvi  ’  and 
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66  is  the  angle  of  wing  rotation. 

Expression  (17)  represents  an  upper  bound  on  ihe  actual  value  of  Fy/i 
(employing  the  upper -bound  principle  discussed  on  page  3  ).  To  obtain 
^  lo we r  bound ,  consideration  is  given  to  the  equilibrium  of  the  joints 
in  the  airframe  truss.  In  satisfying  equilibrium  at  each  joint,  a  re¬ 
duction  in  force  magnitudes  may  be  required  so  as  not  to  exceed 
critical  leads  in  either  tension  or  compression  in  any  one  member. 
Moving  from  joint  to  joint,  a  new  set  of  is  found  with  values  that 
are  less  than  or  equal  to  the  critical  values  and  that  also  satisfy  joint 
equilibrium.  With  these  values,  again  using  expression  (17),  a  lower 
bound  on  Fwi  is  found.  Since  the  lower-  and  upper -bound  values  of 
Fwi  are  not  likely  to  differ  significantly,  their  average  is  then  taken 
as  a  reasonable  value.  See  the  Appendix,  pages  30,  31,  and  32,  for 
computations  on  the  nose-and-wing  impact  configuration. 

GEAR -ONLY  CONTACT.  WITH  VERTICAL  VELOCITY  OTsiT.Y 
CONSIDERED  - - 

A  schematic  diagram  for  this  impact  configuration  is  shown  below. 


FUSELAGE 


OCCUPANT  AND 
SEAT-CUSHION 
COMBINATION 


The  main  gear  is  treated  as  a  collapsible  spring  (in  general,  non¬ 
linear).  The  airframe  is  considered  to  remain  essentially  intact  for 
this  type  of  impact,  as  serious  spinal  injury  could  resul:  from  a 
vertical  deceleration  considerably  less  than  that  required  to  collapse 
the  cabin.  A  standard  seat-cushion  combination  is  assumed,  which 
is  simulated  by  a  linear  spring  with  a  spring  constant  kg  (force  per 
deflection). 

Considering  the  occupant-seat  system,  the  following  energy  balance 
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,11 


may  be  written  for  the  impact  of  the  fuselage  with  the  ground: 


Loss  in  kinetic  energy  of  man  =  increase  in  elastic  energy  of 

seat  and  cushion, 

or 


where 

mni  mass  of  occupant 

Vi  velocity  of  fuselage  upon  its  contact  with  the  ground 

A  total  deflection  of  seat -cushion  combination  associated 

with  a  critical  vertical  deceleration  of  the  occupant. 


Therefore, 


(18) 


To  obtain  the  original  vertical  velocity  of  the  aircraft  prior  to  contact 
of  the  gv;ar  with  the  ground,  an  energy  balance  associated  with  the  gear 
collapse  phase  may  be  written  as  follows: 

2  *  2 

where 

Vq  is  initial  vertical  velocity  of  aircraft 

Ug  is  energy  required  to  collapse  (completely)  the  main  gear. 
Hence, 


(19) 
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Thus,  assm.'ing  that  the  force -displacement  relation  for  the  main  gear 
is  known,  Ug  naay  be  determined  (the  area  under  the  force -displacement 
curve).  Then,  substituting  for  Vj  from  equation  (18),  the  critical  ver¬ 
tical  velocity  Vq  for  this  impact  configuration  is  calculated  by  equation 
(19).  See  the  Appendix,  pages  33,  34,  and  35  for  sample  computations 
on  the  gear -only  type  of  impact  configuration. 
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DISCUSSION 


The  objective  of  the  present  analysis,  as  stated  earlier,  is  to  provide 
a  basis  for  evaluating  proposed  structural  modifications  from  the 
crashworthiness  viewpoint  or  to  indicate  design  possibilities  that  would 
improve  crashworthiness.  Changes  in  structure  might  take  the  form 
of  any  of  the  following: 

1.  General  strengthening  of  a  given  bay  to  increase  its  collapse 
load. 

2.  Increase  in  energy-absorption  capacity  of  a  given  bay  without 
increasing  the  collapse  load. 

3.  General  geometric  changes,  such  as  lengthening  a  bay  or  re¬ 
locating  the  cabin. 

4.  Change  in  component  design,  such  as  a  change  in  the  type  of 
the  main  landing  g<'ar. 

5.  Seat  or  seat  tie-down  changes. 

6.  Minor  structural  changes,  such  as  strengthening  a  single  tube 
or  adding  a  brace. 

7.  Changes  designed  to  reduce  effective  mass  of  the  aircraft 
upon  impact. 

8.  Cabin  modifications  to  increase  cabin  strength  or  perhaps  to 
produce  a  more  favorable  mode  of  collapse. 

I'ach  of  these  is  considered  separately  below,  in  the  light  of  the  analysis 
and  computations  made  so  far  (reference  the  Appendix). 

1.  General  Strengthening  of  a  Given  Bay 

For  a  nose -only- contact  impact  configuration,  an  increase  in  the  first- 
bay  collapse  load  would  increase  the  critical  impact  velocity,  thus 
reducing  the  probability  of  cabin  collapse  during  an  accident.  For 
example,  if  Pj  (for  the  Piper  J-3  investigated)  were  increased  by  50 
percent,  the  critical  velocity  would  increase  from  34  miles  per  hour 
to  45  miles  per  hour  provided  the  first-bay  increase  in  strength  were 
accomplished  without  overloading  the  individual  structural  members 
of  the  cabin.  However,  such  a  change  would  also  increase  the  decel- 
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fe_££_Uon  experienced  by  the  cabin  floor  from  30.  2  G  to  approximately 
45  G.  Hence,  depending  up^^n  the  relative  significance  of  the  two  effects 
produced  (one  favorable  and  one  adverse),  the  increase  in  Pj  may  or 
may  not  be  a  change  that  improves  crashworthiness  significantly. 

2.  Increase  in  Energy  Absorption  Capacity  of  a  Given  Bay 

Without  changing  the  force  required  to  initiate  collapse  of  the  first 
bay,  the  energy  absorbed  by  the  first  bay  may  be  increased  by  im¬ 
proving  the  shape  of  the  load -deflection  curve.  For  example,  if, 
through  additional  structure  or  the  use  of  foamed  plrstic  reinforcement, 
the  average  collapse  load  were  increased  by  50  percent,  the  critical 
velocity  for  the  nose -only  impact  (of  the  Piper  J-3)  would  increase 
from  34  miles  per  hour  to  45  miles  per  hour.  Since  the  maximum 
force  on  the  first  bay  remains  unchanged,  the  deceleration  experienced 
by  the  cabin  floor  would  remain  the  same. 

3.  General  Geometric  Changes 

If  the  first  bay  were  lengthened  by  50  percent  while  the  maximum 
collapse  force  was  maintained,  the  critical  velocity  for  the  nose-only 
impact  (Piper  J-3)  would  increase  from  34  miles  per  hour  to  39  miles 
per  hour.  No  change  in  cabin  deceleration  would  accompany  this. 

4.  Change  in  Component  Design 

If  a  change  were  made  in  t.he  main  landing  gear  so  as  to  increase  its 
energy -absorbing  capacity  (for  collapse)  by  a  factor  of  two,  then  the 
critical  vertical  velocity  (Piper  J-3)  would  be  increased  from  10  miles 
per  hour  to  12.2  miles  per  hour. 

5.  Seat  an  1  Seat  Tie-Down  Changes 

The  present  investigation  is  not  concerned  with  the  relationship  between 
occupant  tie-down  effectiveness  and  survivability;  however,  in  the 
vertical  type  of  impact  (main  gear  only),  the  stiffness  of  the  seat  and 
cushion  combination  enters  into  the  analysis.  The  relative  improve¬ 
ment  in  crashworthiness  effected  by  minor  seat  changes  as  against 
major  structure  changes  is  of  interest.  If,  for  example,  the  seat  and 
cushion  st.'ffness  were  reduced  by  a  factor  of  two,  the  critical  vertical 
velocity  would  increase  from  10  miles  per  hour  to  12.2  miles  per  hour. 
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6.  Minor  Structural  Changes 


It  does  not  appear  that  a  minor  structural  cliange  such  as  strengthening 
a  single  tube  or  adding  a  brace  will  significantly  improve  crashworthi¬ 
ness.  The  possible  exception  lies  in  the  cabin  structure,  which  is 
treated  later. 

The  tubular  structure  airframe  is  generally  redundant  to  the  point  that 
the  strength  of  any  one  tube  v/ill  not  greatly  affect  the  collapse  load  of 
a  bay  or  section  of  the  airframe. 

7.  Changes  Designed  to  Reduce  Effective  Mass  of  Aircraft  upon  Impact 

If  the  effective  mass  of  an  aircraft  were  reduced  upon  impact  by  design 
changes  (such  as  design  features  to  cause  wings  and  tail  section  :o 
separate  from  the  cabin),  a  reduction  in  critical  velocity  could  be  ef¬ 
fected  (for  nose -only-impact  configuration).  For  example,  if  it  were 
possible  to  reduce  the  effective  mass  by  one-third  (in  the  Piper  J-3 
during  nose  impact,  the  critical  velocity  would  increase  from  34  miles 
per  hour  to  42  miles  per  hour.  At  the  same  time,  the  cabin  decelera¬ 
tion  would  increase,  from  30.  2  G  to  45  G.  Again,  such  a  change  pro¬ 
duces  both  beneficial  and  adverse  effects,  and  the  net  benefit  attained 
would  depend  upon  the  relative  significance  of  these  effects. 

8.  Cabin  Modification 

Although  not  investigated  in  detail  during  this  analysis,  improved 
crashworthiness  appears  to  be  probable  from  changes  in  cabin  struc¬ 
ture.  If  the  cabin  (defined  as  the  protective  5hel)  immediately  sur¬ 
rounding  the  occupants)  were  made  n-  :ch  stronger  than  other  airframe 
structure,  impact  forces  could  be  transmitted  through  the  cabin  and  the 
energy  absorbed  in  deformation  elsewhere.  The  present  analysis  has 
considered  cabin  strength  nly  to  verify  that,  for  the  particular  forces 
associated  with  a  given  collapse  mode,  the  cabin  would  remain 
essentially  intact. 
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LIMITATIONS  OF  ANALYSIS 


The  analysis  made  so  far  contains  both  limitations  and  uncertainties. 
As  it  was  developed  for  the  tube  and  fabric  type  of  light  aircraft  and 
was  in  particular  oriented  toward  the  Piper  J-3  type  of  construction, 
it  would  be  inadequate  for  any  aircraft  that  departs  greatly  from  the 
Piper  J-3  in  structure.  Moreover,  the  analysis  contains  several 
arbitrary  assumptions  that  deserve  further  investigation.  Static 
loading  is  assumed  on  any  given  substructure.  An  intermediate  value 
of  column-end  fixity  is  assumed.  Simplifying  assumptions  were  made 
in  connection  with  energy  dissipation.  In  the  computations  (reference 
the  Appendix),  values  were  assumed  both  for  the  soil  elastic  constant 
and  for  the  effective  mass  of  the  aircraft.  These  assumptions  pre¬ 
sently  stand  as  vmcertainties . 
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RESULTS  OF  NUMERICAL  EXAMPLES 


The  following  table  summarizes  the  results  of  computations  (reference 
the  Appendix)  for  the  three  impact  configurations  discussed  using  the 
Piper  J-3  as  an  example  aircraft  design.  The  critical  velocity  is 
understood  to  mean  that  velocity  (for  the  mathematically  idealized 
conditions)  that  produces  airframe  collapse  just  short  of  cabin  collapse, 
or,  in  the  case  of  vertical  velocity,  that  produces  a  vertical  decel¬ 
eration  of  20G. 


Impact 

Configuration 

Critical 
Velocity  (mph) 

Maximum 
Deceleration  (G) 

Nose -only  contact 

34.  0 

30.2 

Nose-and-wing  contact 

35.7 

33.6 

Gear-only  contact 

10.0 

20.0 

(vertical  velocity) 


APPENDIX.  NUMERICAL  EXAMPLES 


The  computations  for  the  Piper  J-3,  associated  with  equations  3,  4,  7, 
8,  9i  iO,  13,  15,  16,  17,  18,  and  19,  appear  on  the  attached  data 
sheets. 

Tube  data  were  taken  from  Piper  structural  drawings.  The  effective 
mass  for  the  entire  aircraft  was  arbitrarily  taken  as  90  percent  of  the 
maximum  gross  weight.  This  is  based  upon  presumed  flailing  of 
occupants,  shifting  of  luggage,  and  sloshing  of  fuel.  The  soil  spring 
constant  (3000-pounds -per-inch  penetration)  is  based  upon  reported 
penetrations  and  force-niass -acceleration  data  of  the  NACA  crash 
tests  using  Piper  J-3  aircraft  which  were  crashed  into  tamped  earth 
embankments.  The  gear  collapse  energy  for  the  Piper  J-3  may  be 
computed  from  elastic  data  on  stretch  of  bungee  cords  (four  loops). 

Ihe  stretch  as  installed  on  the  aircraft  is  4  inches,  and  collapse 
stretch  is  an  additional  25.  2  inches. 
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FIRST  BAY  FORCE  COMPUTATIONS 


Piper  J-3 


Member 

^i 

yi 

*i 

?i 

*i 

(Jc  -  X 

(9  -  y 

(2  -  z 

0 

17.0 

17.0 

1 

9.75 

13.5 

3.75 

-1.5 

-25.0 

-23.5 

0 

26.2 

2G.2 

2 

9.75 

12.25 

2.5 

-1.5 

-25.0 

-23.5 

0 

26.2 

26.2 

3 

9.25 

0 

-9.25 

-20.0 

-25.0 

-5.0 

0 

26.2 

26.2 

4 

9.75 

12.25 

2.5 

-20.0 

-25.0 

-5.0 

0 

17.0 

17.0 

5 

0 

13.5 

13.5 

-20.0 

0 

20.0 

Fy  cri  1 


305 

17.5 


1245 

35.3 


797 

28.2 


718 

26.8 


871 

29.5 


.971 

.085 


.927 

-.177 


.977 

-.186 


.576 

.677 


.59 

.035 

16.5 

.715 

.035 

33.0 

.59 

.035 

26.0 

.59 

.035 

25.0 

.653 

.035 

26.0 


6.12 

4.54 


2.73 

5.50 


2.47 

4.54 


2.67 

4.54 


3.34 

5.02 


4.41 


2.03 


2.29 


2.61 


*For  the  purpose  of  coinputing  buckling  loads,  the  effective  length 
of  the  tube,  Le,  was  used,  th'’  length  of  uniform  cross-section  of 
the  tube. 


7  I  =  13.26 


=  26.5  kips 

Maximum  deceleration  for  nose-only  contact  (assuming  effective 
mass  of  .9  of  actual  mass): 


Pi  =  m’a 


m' 

a  = 

*  7173  ^ 

a  =  30.4  G 


m'  = 


(.9)(.97)  .  . 
- g - kips 
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Use: 


K  “  3  kips/ inch  (based  upon  data  from  NACA  crash  tests) 
s  =  17  inches 
m'  =  (.9)(970)  =  873  lb. 
m  =  (.9)(1200)  =  1080  lb. 


Therefore, 


Z2 

Z 

v2 


V 


/  Pi  m  \  , 

V  + 1 

[l  ^K?/  ' 

[  2  +  .52  (1.235) 

1  Z  +  (1  ♦ 

1.32  +  ^1.74  -  1 

.52  =  1.79 

(26.5)(1.79)(17)(32.2) 

"in^  ■  (.8>3j(12) 

2478 

49,8  fps  ^  34  nph 
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•Syo  =  *yr  = 


=  fSB 


Therefore,  6B  = 


=  2b6B 


^  2b(6B  +  66) 


=  2b(l  +  t)  66 


An  upper  bound  to  the  external  load,  required  to  sustain  the 

assumed  mechanism  may  be  found  by  en^loying  the  virtual  work  principle. 

6 (External  work)  *  6 (Internal  work) 


F„,r66  =  T  T  6. 

wi  *•  crj^  1 


F„ir66  =  F^,d^(6xB  ♦  6xc) 


=  '■cdjf-**®  ♦  2b(i  ♦l)«e] 

^  ^  2b  [2a  ,1 

f^wl  =  ^cd„  r  If  * 


=  (2.79)  (2.7)  =  1.04  kips 


The  upper  bound  would  also  be  a  lower  bound  if  equilibrium  of  internal 
forces  were  satisfied  at  each  joint.  A  consideration  of  each  joint 
free  body  shows  that  equilibrium  can  be  satisfied  by  appropriate 
internal  forces  compatible  with  the  load,  F^^j,  found  above.  Thus,  F„j 
is  correct. 


I 


-  [  2  +  .666  -  .020  ]  Y  +  1.52  =  0 
=  1.322  yT.75  -  1.52  =  1.80 


Therefore, 


^  ^  ^wi)sY 


1.64 

my  =  •—  slugs 


v2  = 


(2.74)(17)(1.80)(32.2) 

“TTHTTOTO^t^ — 


=  2749 


V  -  52.4  fps  or  35.7  mph  . 


Deceleration  for  Nose  and  Wing  Contact 

__  n  _  ^wl  _  26.5  t  1.04  _  ,  . 

Equation  (13)  F„  i.84  kips 

«w 


Haxivum  possible  force  = 


^wi  ^w2 


=  26.5  t  1.04  +  1.84  =  2.94  kips 
=  J-  .  29^  -  33  7G 

bT  7573 
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Main-Landing-Gear-Contact  Configuration 

Vfhere 

Wj„  is  the  occupant  weight  =  170  lb, 
kg  is  the  assumed  seat  elastic  constant  =  1700  psi 
For  maximum  deceleration  of  20g, 


F  =  ma 

kgA  =  (20g) 

'O' 

or  t  -  2  inches. 

Then,  using  Equation  (18), 


A  /  (1700)(12)(32.2)  7^ 

’  V  170  12 

Vj  s  10,35  fps  , 

For  the  Piper  J-3,  the  energy  of  main  gear  collapse  may  be  obtained 
from  the  load  deflection  curve  of  the  bungee  cord  used: 


PqX  + 


kx^ 

2 


times  number  of  bungee  cords  stretched 


energy  per  cord 


where 


3 


Pq  is  initial  tension  in  cord  =  100  lb, 

X  is  total  elongation  of  each  cord  =  25,2  in. 


k  is  cord  spring  constant  (from  load  deflect icw:  curve)  = 
9,16  lb-in,  elongation. 

(Each  cord  loop  is  28  inches  long  unstretched). 
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Therefore,  for  the  Piper  J-3, 

U  -[(J-00)(35.2)  ^  (9.16)(25.2)^ 
g  I  12  24 

Ug  =  1812  ft-lb. 


Vq  =  14.65  fps 
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Load-Deflection  Cuive  for  Piper  J-3  Landing  Gear  Bungee  Cord 
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